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ABSTRACT 



Pl,: 

6 

We present measurements of [Fe/H] for six M supergiant stars and three giant stars 
within 0.5 pc of the Galactic Center (GC) and one M supergiant star within 30 pc of 



the GC. The results are based on high-resolution (A/ A A = 40,000) K— band spectra, 
taken with CSHELL at the NASA Infrared Telescope Facility. We determine the iron 
abundance by detailed abundance analysis, performed with the spectral synthesis pro- 
gram MOOG. The mean [Fe/H] of the GC stars is determined to be near solar, [Fe/H] 
= +0.12 ± 0.22. Our analysis is a differential analysis, as we have observed and applied 
the same analysis technique to eleven cool, luminous stars in the solar neighborhood 
with similar temperatures and luminosities as the GC stars. The mean [Fe/H] of the 
solar neighborhood comparison stars, [Fe/H] = +0.03 ± 0.16, is similar to that of the 
GC stars. The width of the GC [Fe/H] distribution is found to be narrower than the 
width of the [Fe/H] distribution of Baade's Window in the bulge but consistent with the 
width of the [Fe/H] distribution of giant and supergiant stars in the solar neighborhood. 
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Subject headings: Galaxy: abundances, Galaxy: center, stars: abundances, stars: late- 
type 

1. Introduction 

Morris & Serabyn (1996) and Serabyn & Morris (1996) have summarized the relative impor- 
tance of star formation, gas inflow, and gas outflow in the Galactic Center (GC). Gas from the 
inner disk flows into the nucleus, perhaps driven by a bar potential (Stark et al. 1991; Morris & 
Serabyn 1996), and some of the stellar mass-loss from the bulge may fall into the GC (Blitz et al. 
1993; Jenkins & Binney 1994). Serabyn & Morris (1996) have proposed that this inflow of disk gas 
results in sustained star formation in a "central molecular zone" (CMZ) within a radius from the 
GC of R < 200 pc. The GC is observed to be currently forming stars in the CMZ at a rate of 
0.3-0.6 Mq yr _1 (Glisten 1989). Some of the GC gas incorporated into stars, enriched by stellar 
nucleosynthesis, will be returned to the GC interstellar medium by stellar mass-loss and super- 
novae. X-ray observations of hot gas in the GC suggest that some gas could be driven temporarily 
or permanently from the GC by a galactic fountain or wind (Blitz et al. 1993; Morris & Serabyn 
1996), however, making it unclear how much enriched material is incorporated into following gen- 
erations of star formation. The presence of strong magnetic fields in the GC also likely plays a role 
in whether enriched gas is driven from or is retained within the GC. 

Continued star formation in the central few hundred parsecs of the Galaxy may lead to higher 
metallicities within the CMZ. Chemical abundances in the disks of spiral galaxies are observed to 
reach their highest values at the center (see Pagel &: Edmunds 1981; Shields 1990). H II regions, 
planetary nebulae, and OB associations in the Milky Way also show a radial metallicity gradient 
(Shaver et al. 1983; Fich & Silkey 1991; Maciel & Koppen 1994; Simpson et al. 1995; Vflchez & 
Esteban 1996; Rudolph et al. 1997; Afflerbach et al. 1997; Smartt & Rolleston 1997; Gummersbach 
et al. 1998). Chemical evolution models strive to explain these gradients by considering the relative 
star formation and gas infall/outflow rates, and the metal abundance of the gas compared to the 
stars (see Audouze & Tinsley 1976; Rana 1991). 

The metallicity of GC stars is needed in order to constrain models and understand the stellar 
processes in the central parsecs of our Galaxy. Measurements of stellar metallicities in the center of 
the Milky Way Galaxy, which is obscured by approximately 30 mag of extinction at V, are only now 
beginning to be feasible through infrared studies. Carr et al. (2000) found a solar Fe abundance 
in IRS 7, an M supergiant at a distance of R = 0.2 pc from the GC, from a detailed abundance 
analysis of CSHELL (A/AA = 40,000) K— band spectra. Our work analyzes similar data for several 
M supergiant and giant stars located in the central cluster, within 0.5 pc of the GC, and also one 
M supergiant star located in the Quintuplet cluster, 30 pc away from the GC. 
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2. Observations 

2.1. Sample Selection 

The GC stars were selected with three requirements. First, they have to be brighter than K 
= 9.5, a limit set by the sensitivity of the instrument. An hour of integration time is required 
for high signal-to-noise (S/N) spectra for a K = 9 star when observing with CSHELL at the 
NASA Infrared Telescope Facility (IRTF). Second, the stars should be cool luminous stars, to 
ensure we are observing photospheric absorption lines rather than the stellar wind emission lines 
characteristic of hot luminous stars in the GC. Third, these cool luminous stars should have little 
or no water absorption. This is because published line lists for water (wavelengths, excitation 
potentials, dissociation potentials, oscillator strengths, and damping constants) are not yet adequate 
for detailed abundance analysis of high-resolution stellar spectra. The amount of water absorption 
is determined from H and K— band low-resolution spectra (A/AA ~ 500) taken with the IRS at 
CTIO (Ramirez et al. 1999). The selected sample consists of nine stars located in the central cluster 
(R < 0.5 pc), and one star located in the Quintuplet cluster (R = 30 pc). The central cluster stars 
are IRS 7, IRS 11, IRS 19, IRS 22, BSD 72, BSD 114, BSD 124, BSD 129, and BSD 140 (names 
and positions are given by Blum et al. 1996a). The star located in the Quintuplet cluster is VR5-7 
(Moneti et al. 1994). Eleven cool, luminous stars in the solar neighborhood were selected from 
the literature to be observed and analyzed the same way as the GC stars. These eleven stars 
have known abundances from detailed analysis in the optical (Luck 1982a,b; Lambert et al. 1984; 
Smith & Lambert 1985, 1986; Luck & Bond 1989) and in the infrared (Carr et al. 2000). These 
cool, luminous stars in the solar neighborhood were carefully selected to be in the same range of 
effective temperature and surface gravity as the GC stars. This is very important in order to allow 
a differential analysis comparison, which will cancel possible systematic errors such as NLTE effects 
(see Sec. 4.2). The selected cool, luminous stars in the solar neighborhood are listed in Table 1. 



2.2. Data Acquisition and Reduction 

The K— band high-resolution (A/AA = 40,000) spectra used in the abundance analysis were 
taken with CSHELL at the IRTF. The IRTF is a 3-m telescope located at Mauna Kea Observatory 
in Hawaii, and CSHELL is the facility cryogenic infrared echelle spectrograph (Tokunaga et al. 
1990). The observations were carried out in 1993 May, 1996 June, 1996 August, 1997 July, 1998 
May, 1998 July, and 1999 June. The detector for the 1993 observations was a 256x256 NICMOS-3 
HgCdTe array, which provided a spectral coverage of ps 1000 km s _1 or 73 A. The detector for the 
rest of the observations was a 256x256 SBRC InSb array, which gives a smaller coverage of « 750 
km s _1 or 55 A. The slit used to achieve a resolution of A/AA = 40,000 is 0.5" wide, 30" long, 
with a pixel scale of 0.25" per pixel for 1993 data and 0.20" per pixel for later data. Our highest 
quality data were taken in 1998 and in 1999 when new order-separating circular variable filters were 
installed, providing spectra free of fringes. 
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Eight iron lines were selected by inspecting the high resolution atlas of cool stars (Wallace 
& Hinkle 1996) and by analyzing synthetic spectra, which includes atomic and molecular features 
(see details in Section 3.0 and in Carr et al. (2000)). Thousands of lines of molecules like CN, 
CO and H2O are present throughout the infrared spectrum of cool stars. The selected iron lines 
are almost free of known molecular contamination and unidentified lines. The eight iron lines are 
observed using three grating settings of CSHELL. Atomic parameters of the eight iron lines used 
in the abundance analysis are listed in Table 2. 

Spectra were acquired at two different positions along the slit and reduced separately, to 
determine any systematic effects like fringes. Several individual spectrum pairs were observed per 
star, until the desired signal to noise ratio was reached. Stars of spectral type A or B were observed 
as close to the program star's airmass as possible to correct for telluric absorption features and to 
remove fringes. Such stars have no significant spectral features in the observed wavelength regions. 

Image Reduction and Analysis Facility (IRAF) and VISTA were used for data reduction. The 
reduction process began by flat fielding the individual spectra with calibration lamp flats, which 
were taken for every grating setting. Sky subtraction was performed by subtracting two consecutive 
spectra taken at different positions along the slit. Bad pixels were replaced by an interpolated value 
computed from neighboring pixels in both the dispersion and perpendicular directions. Individual 
spectra were extracted with a 5 pixel wide aperture using the APSUM package in IRAF. Extracted 
spectra observed at the same position along the slit were averaged to produce two spectra per star, 
one spectrum for each position along the slit. Spectra of the program stars were then divided by 
the appropriate A or B type atmospheric standards, observed at the same position along the slit 
and reduced in the same way, to remove telluric absorption features and fringes. Each spectrum 
was then multiplied by a blackbody of the same temperature as the A or B star to put the spectra 
on a relative flux density (F^) scale. The temperature of the A or B star was determined from its 
spectral type. More than five telluric absorption lines were used to obtain wavelength solutions. The 
wavelength calibration was performed using the tasks of VISTA. Spectra were shifted in wavelength 
to correct for radial velocity differences, by comparing the observed and synthetic spectra. The 
spectra were normalized by a linear fit of continuum bands. These continuum bands were chosen 
by comparing the observed spectrum of a Ori from the high-resolution atlas (Wallace & Hinkle 
1996) with the synthetic spectrum of a Ori. The bands selected are regions in the spectrum where 
the synthesis and the normalized observations of a Ori are coincident and equal to unity. We 
estimated the error per pixel as the difference of observations taken at two positions along the slit. 
We computed the mean signal to noise ratio (S/N) per pixel by the average of the signal to noise 
ratio per pixel for all spectra observed for that star. The mean S/N per pixel for each star is listed 
in Tables 1 and 3. 
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3. Abundance Analysis 

The abundance analysis was done using a current version of the LTE spectral synthesis program 
MOOG (Sneden 1973). The program requires a line list giving the wavelength, excitation potential, 
gf-values, and damping constants for all atomic and molecular lines that contribute to the spectrum. 
The method used to construct the line list is described in Sec. 3.1. In addition, an input model 
atmosphere for the effective temperature and surface gravity appropriate for each star and a value 
for the microturbulent velocity is also required. The solar-abundance model atmospheres from Plez 
(Plez 1992; Plez et al. 1992) were used in our analysis because they are the only ones that include 
gravities and temperatures as low as those of the Galactic Center stars. The Plez models also 
include sphericity, which is appropriate for supergiant stars (Plez et al. 1992). MOOG, however, 
does not account for sphericity in its calculations. Carr et al. (2000), who also used MOOG, 
compute the effective temperature and microturbulent velocity of a Ori, VV Cep (HR8383), and 
(3 And using different sets of model atmospheres, including the Plez models and two plane-parallel 
atmospheres. When the stellar parameters derived from each model were used self-consistenly to 
derive [Fe/H], the result for [Fe/H] was the same. The determination of the stellar parameters for 
our sample is discussed in Sec. 3.2. 

3.1. Atomic and Molecular Parameters 

The line list was created the same way as Carr et al. (2000), using two main steps: 1. The solar 
spectrum was used to determine the gf-values and the damping constants. 2. Minor adjustments 
were made by comparison to the Arcturus spectrum (Wallace & Hinkle 1996). An initial line list 
was compiled using Fe I line positions and energy levels from Nave et al. (1994) and the Kurucz 
(1993) line lists for CN and other atomic lines. The damping constants for all atomic lines were 
initially set to twice that of the Unsold approximation for van der Waals broadening (Holweger et al. 
1991). A synthetic spectrum was generated for the Sun, using the Kurucz (1993) solar model. The 
gf-values and damping constants were adjusted to match the observed solar spectrum (Livingstone 
Sz Wallace 1991), when needed, and unidentified lines were noted. Then, a synthetic spectrum 
was generated for Arcturus, using the model atmosphere, stellar parameters, and abundances from 
Peterson et al. (1993). The solar line list provided a good match to the Arcturus spectrum. The 
gf-values of low excitation lines, not observable in the Sun, had to be adjusted. The final Fe I 
atomic parameters are listed in Table 2. A synthetic spectrum of a Ori, generated using the final 
line list, is compared in Figure 1 to the observed spectrum from Wallace & Hinkle (1996) in the 
region of each of our eight Fe I lines. A Plez model atmosphere with the stellar parameters listed 
in Table 3, and the CNO abundances from Lambert et al. (1984) were used. 

The depth of line formation was examined using contribution functions, which provide an 
indication of spectral line formation coming from different layers of a stellar model atmosphere (see 
details in Edmonds 1969; Sneden 1973). The determination of the depth of line formation is needed 
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to make sure each line is formed in the range of opacity (r) covered by our model atmosphere. The 
eight iron lines used in this abundance analysis are formed in the range of -2.0 < log r < 0.3, where 
Plez model atmospheres cover a range of -5.6 < log r < 2.6 . Note that the Sc I lines in supergiant 
stars, which neither we nor Carr et al. (2000, who included hyperfine splitting calculations for Sc I 
lines) are able to model correctly, are shown by the contribution function to be partially formed at 
optical depths outside the range of the Plez models. We therefore approximately model Sc I lines 
by scaling the Fe I line profiles to the correct depth and width to allow us to separate the one Fe I 
line at 22399 A which is blended with a Sc I line at 22400 A. In the case of IRS 7 this procedure 
was not successful. 



3.2. Stellar Parameters 

3.2.1. Effective Temperature 

The effective temperature (T e g) is a key parameter for any abundance analysis. For the 
cool, luminous stars in the solar neighborhood, the effective temperatures were taken from the 
literature. Some effective temperatures come from spectroscopic measurements (Luck 1982a, b; 
Luck & Bond 1989; Carr et al. 2000). Other effective temperatures come from the calibration of 
(V — K) colors vs. T e g, where T e g- comes from angular diameters measured by lunar occultation 
(Smith & Lambert 1985, 1986). a Ori has an T e g determination from its angular diameter measured 
by lunar occultation (T e ff = 3605 ± 43 K, Dyck et al. 1998), which is very good agreement with 
the spectroscopic value obtained by Carr et al. (2000) (T e g- = 3540 ± 260, see Table 1). 

For IRS 7, the brightest infrared source in the GC, Carr et al. (2000) determined T e g spectro- 
scopically by requiring the carbon abundance derived from CO lines to be independent of the lower 
excitation potential. We adopt the stellar parameters for IRS 7 from Carr et al. (2000), which were 
derived using the standard Plez grid of model atmospheres. Their T e g analysis relied largely on 
the second overtone CO lines present in the H— band (1.5/um - 1.8^m), because most of the first 
overtone CO lines present in the K— band are saturated. 

For the remaining nine stars in the GC the H— band spectra are too faint to be observed 
with CSHELL at the IRTF, because of high extinction. Therefore CO cannot be used to compute 
the effective temperature. For those nine GC stars, broad molecular bands present in H and 
K— band low resolution spectra (A/AA ~ 500) are used to estimate T e g- (Ramirez et al. 1999). 
These features are CO (2.3 //m) and H2O (1.9 //m). The CO strength increases with decreasing 
T e g and decreasing surface gravity (log g). The H2O strength also increases with decreasing T e ff , but 
increases with increasing log g. These two features together, therefore, provide two-dimensional 
spectral classification (Baldwin et al. 1973; Kleinmann & Hall 1986; Blum et al. 1996b). The 
low resolution spectra of GC stars were taken with the IRS at CTIO, used in cross-dispersed 
mode to acquire the H and K— band simultaneously (Ramirez et al. 1999). Nearby late- type stars 
with T e g determinations based directly or indirectly on the lunar occultation technique (Smith & 
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Lambert 1985, 1986, 1990; Fernandez- Villacahas et al. 1990; McWilliam 1990; Dyck et al. 1998; 
Richichi et al. 1998) were observed to calibrate our T e g- determination. The low resolution spectra 
of cool, luminous stars stars with known T e ff were taken with the IRS at CTIO, and with TIFKAM 
at MDM Observatory (Ramirez et al. 1999). The distinction between giants and supergiants is 
determined by the presence or weakness of H 2 0, respectively. The CO index is computed in 
the same way as in Blum et al. (1996b). Once the distinction between giants and supergiants is 
established, T c g is determined from the relations shown in Figure 2, which are the best unweighted 
linear fit to the data. VR 5-7, IRS 19, IRS 22, BSD 72, BSD 124, and BSD 129 were classified as 
supergiant stars, and IRS 11, BSD 114 and BSD 140 were classified as giant stars. The luminosity 
class for GC stars is also listed in Table 3. Our values for T e g- for GC stars agree with the results 
from Blum et al. (1996b), for the stars that are common to both samples. Our value for T e g- from 
the CO index of IRS 7, 3400 ± 300 K, is also in good agreement with the spectroscopic value 
of T e g from Carr et al. (2000), 3470 ± 250 K. IRS 7 has a low carbon abundance (Carr et al. 
2000, [C/H] = -0.8). The good agreement between the T e g from CO index and the spectroscopic 
T e g suggests that our derived values of T e g are not strongly affected by moderate variations in the 
abundance of carbon. The typical uncertainty in T e g- computed from the CO index is 300 K for 
supergiants and 280 K for giants. The T e R for cool, luminous stars in the solar neighborhood is 
listed in Table 1 and for GC stars is listed in Table 3. 



3.2.2. Surface Gravity 

For the cool, luminous stars in the solar neighborhood, the surface gravity, g, from the literature 
was used (Luck 1982a,b; Smith & Lambert 1985, 1986; Luck & Bond 1989; Carr et al. 2000). For 
the GC stars, the surface gravity is determined from the relation: 

log g = log(M/M ) + 4 log(T cff /T cff0 ) - log(L/L Q ) - log g e . (1) 

T e ff has been already determined. The luminosity, L, is computed from Mbd, which is determined 
from K, Ak, the GC distance, and the bolometric correction, B.C. A GC distance of 8 kpc is 
assumed (Reid 1993), K and Ak are from Blum et al. (1996a), and the B.C. is taken from Elias 
et al. (1985). The final Mbd for the GC stars is listed in Table 3. The uncertainty in Mbd from 
the GC stars comes mainly from uncertainties in the extinction curve, E(H — K)/Ak , and it 
estimated to be ±0.4 (Carr et al. 2000). Once Mboi and T e ff are determined, the stars are placed 
in the HR diagram, and masses, M, are obtained by over-plotting evolutionary tracks of varying 
initial masses. Solar metallicity evolutionary tracks (Schaller et al. 1992) are used to determine 
masses for supergiant stars. For the coolest stars, that fall outside the T e g- range of Schaller et al. 
(1992) evolutionary tracks, asymptotic giant branch (AGB) tracks from Marigo et al. (1996) were 
used to estimate their masses. Note that the stellar evolution models assume mass loss for the 
more massive stars, and the current mass of the star at its present age is assumed. Twice solar 
metallicity evolutionary tracks (Schaerer et al. 1993) gave the same result for the derived masses. 
The value of log g for GC stars is listed in Table 3. The uncertainty is estimated considering the 
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uncertainties in mass, effective temperature and luminosity. Figure 3 shows the HR diagram of 
GC stars and also the cool, luminous stars in the solar neighborhood. It is seen that the GC stars 
and the cool, luminous stars in the solar neighborhood occupy the same place in the HR diagram, 
hence are very similar types of stars. 

3.2.3. Microturbulent Velocity 

The values for the microturbulent velocity (£) for the cool, luminous stars found in the literature 
show considerable variation among authors. Carr et al. (2000) found for VV Cep (HR 8383) a much 
lower value (£ = 3.7 ± 0.2 km s" 1 ) than the one from Luck (1982b, £ = 5.0 ± 0.5 km s" 1 ). Carr et 
al. (2000) used CO lines of the same excitation potential, for which the carbon abundance should 
be independent of the equivalent widths of the lines when £ is correct. Luck (1982b) used the same 
principle, but used Fe I lines instead of CO lines. 

Seven of our Fe I lines have similar excitation potential, and those lines were used to derive 
the microturbulence for the cool, luminous stars in the solar neighborhood. This was done in an 
iterative way. An Fe abundance was derived from the synthetic model of cool, luminous stars in 
the solar neighborhood, which included lines of Fe, CN, and other lines. Using the derived Fe 
abundance, a synthetic model was computed containing only Fe lines, and an equivalent width 
was measured for the Fe lines. This way the equivalent width is free from contributions from CN 
or other species. Then, the derived Fe abundances from the synthesis were plotted vs. the Fe 
equivalent widths. This process was repeated with different microturbulent velocities until the Fe 
abundance was independent of the Fe equivalent width. The final equivalent widths of the Fe lines 
used in this process are listed in Table 4. The obtained value of £ for cool, luminous stars in the 
solar neighborhood is listed in Table 1. The uncertainty in the obtained microturbulent velocity 
comes from the uncertainty in the slope in the graph of Fe abundance vs. equivalent width, based 
on the scatter of the data points. Our value of £ for a Ori (£ = 2.8 ± 0.2) is slightly different 
from the value of £ from Carr et al. (2000) (£ = 3.23 ± 0.15). There is good agreement between 
our values of £ and the ones given by Smith & Lambert (1985, 1986), where the mean difference is 
0.15 km s _1 . Our values of £ are systematically lower than the values given by Luck (1982b) and 
Luck & Bond (1989) (HR 8383, HD 202380, HD 163428, BD+59 594, HD 232766), where the mean 
difference is 2.4 km s _1 . Our value of £ for HR 8726 is in good agreement with the value given by 
Luck (1982a). 

We do not have sufficient S/N for the GC stars to apply this technique for finding £. Instead, to 
get the microturbulent velocity for the GC stars, a relation between log g and £ is used. McWilliam 
(1990) showed that log g and £ are related in G and K giants. The values of £ and log g for cool, 
luminous stars in the solar neighborhood, including the values from Carr et al. (2000) and Smith & 
Lambert (1985, 1986, 1990), are plotted to get a relation to apply to the GC stars. Figure 4 shows 
this relation for stars with T e g and log g similar to those of the GC stars. An unweighted linear fit 
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to the data gives : 



f = (2.78 - 0.82 x log g) km s 



-i 



(2) 



The uncertainty in the microturbulent velocity is estimated from the uncertainty of the fit 
(±0.4 km s" 1 ) and the uncertainty from the surface gravity. The value of £ for IRS 7 obtained by 
this fit, £ = 3.3 ± 0.4 km s _1 , is in good agreement with the value of £ from Carr et al. (2000), £ 
= 3.30 ± 0.34 km s^ 1 , obtained using CO lines. 



The synthetic spectra have to be convolved with a macroturbulent broadening function (Gray 
1992), and with instrumental broadening, to match the line widths of the observed spectrum. The 
instrumental broadening is a gaussian of full width half maximum given by the spectral resolution 
of the instrument (7 km s" 1 ). CN lines in two bands (21798 A - 21804 A; 22403 A - 22409 A) 
were used to determine the macroturbulent velocity (() for the cool, luminous stars in the solar 
neighborhood. Several synthetic spectra were generated by MOOG using different values of ( 
and N abundances. Then, \ 2 was computed for each synthetic spectrum, and the best set of N 
abundance and £ was chosen for the minimum \ 2 ■ Note that this procedure is not measuring (, 
but determining the best line broadening that fits the data. Also, the N abundance obtained by 
this technique is not real, because the C and O abundances have not been measured independently. 
Figure 5 shows the x 2 contour maps, one for each band, for HR 8726. Figure 5 also shows the 
observed spectrum of HR 8726 and the synthetic spectrum for the values of £ and N abundance 
which give a minimum \ 2 f° r both bands. The macroturbulent velocity for the cool, luminous 
stars in the solar neighborhood is given in Table 1. The uncertainty in £ is given by the standard 
deviation of the values of C derived for each band. This uncertainty is more conservative than the 
uncertainty derived from \ 2 statistics. 

For the GC stars, the same method was used for IRS 7, VR 5-7, IRS 19, IRS 22, and IRS 11. 
For the remaining stars, BSD 72, BSD 114, BSD 124, BSD 129, and BSD 140, only one band with 
CN lines was observed and the band did not have enough S/N to carry out this method. The values 
of £ in Table 1 were used to derive a relationship between log g and £ for supergiant stars. Gray & 
Martin (1979) established a relationship between £ and £ for giant stars. McWilliam (1990) shows 
that £ depends on log g, so by combining these two dependences we derive a relation between ( 
and log g. The values of ( for supergiant stars measured by the broadening of the CN lines are 
plotted vs. log g to get a relation to apply to the fainter GC supergiant stars. Figure 6 shows this 
relation. An unweighted linear fit to the data gives : 



The uncertainty in the microturbulent velocity is estimated from the uncertainty of the fit (±2.4 
km s^ 1 ) and the uncertainty in the surface gravity. Equation (3) was applied to BSD 72, BSD 124, 
and BSD 129. 



3.2.4- Macroturbulent Velocity 



C = (14.5 - 4.40 x log g) km s 



-i 



(3) 
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For giant stars, the range in log g is too narrow to derive a relationship. A mean value of ( = 
9.0 ± 3.0 km s _1 is used for the GC giant stars. This mean value is derived from values of £ for 
solar neighborhood giant stars and IRS 11 measured by the broadening of the CN lines, as given 
in Tables 1 and 3. This mean value was adopted for BSD 114 and BSD 140. 

3.3. Uncertainties 

Our basic observable parameters, T e g- and -Mbd, have uncertainties that come from the tech- 
nique used to compute them (sections 3.2.1 & 3.2.2). The surface gravity uncertainty includes both 
the uncertainty in T e Q and the uncertainty in Mboi (section 3.2.2). The microturbulent velocity 
uncertainty includes both the uncertainties in log g and the uncertainty of the fit used to calculate 
it (section 3.2.3). The same is true for the macroturbulent velocity (section 3.2.4). 

It is important to know what are the effects of each of these uncertainties in the abundance 
determination. The uncertainty in [Fe/H] for each stellar parameter has been estimated by varying 
one of the stellar parameters and computing the difference in iron abundance. A typical uncertainty 
in T e ff of ±300 K results in an uncertainty of =F 0.08 dex in [Fe/H]. A typical uncertainty of ±0.35 
in log g implies an uncertainty of ±0.12 dex in [Fe/H]. A typical uncertainty of ±0.4 km s _1 in £ 
causes an uncertainty of =F 0.09 dex in [Fe/H]. A typical uncertainty of ±1.8 km s _1 in ( translates 
to an uncertainty of ±0.13 dex in [Fe/H]. The iron abundance is not very sensitive to uncertainties 
in T e f[ and £, but it is more sensitive to uncertainties in log g and £. In addition to the uncertainty 
from the stellar parameters, the standard error must also be considered for each star. The standard 
error comes from the scatter in [Fe/H] as derived from individual Fe I lines. All the uncertainties 
from stellar parameters and the standard error for each star are listed in Table 7 and 8. The total 
uncertainty is the quadratic addition of the uncertainties derived from the stellar parameters and 
the standard error in the value of [Fe/H] from individual Fe I lines. For the solar neighborhood stars 
and IRS 7, the quadratic addition of the uncertainties is a good estimation of the total uncertainty, 
because all the uncertainties come from independent measurements. In the cases where £ and £ 
are computed though their relationships with log g, the uncertainties are not independent, and 
the quadratic addition can be an overestimation of the total uncertainty. We estimate that the 
total uncertainty will decrease by no more than 0.05 dex, if we consider the correlation of the 
uncertainties in log g, £ and (. 

Systematic errors may be present as well. They may come from unidentified lines, which could 
make a contribution that is not included in our synthetic spectrum because of the lack of atomic and 
molecular parameters. They may also come from errors in the gf-values determination, or failure 
of the model atmospheres to correctly model the stars we observe, because model atmospheres use 
solar abundance ratios, or NLTE effects. 
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3.4. NLTE effects 

The iron abundance could be affected by departures from LTE. The main NLTE effect in late- 
type stars is caused by overionization of electron donor metals by ultraviolet radiation (Auman & 
Woodrow 1975). Using the Saha equation, we found that in our typical star 93 % of iron is neutral. 
Thus overionization should be smaller in our stars than in warmer giants, where only a few percent 
of iron is neutral. 

It is known from empirical studies that abundances derived from low excitation lines give sys- 
tematically lower abundances (~0.3 dex) than the ones derived from high excitation and ionized 
lines in late-type stars; high excitation and ionized lines give very similar abundance results (Ru- 
land et al. 1980; Steenbock 1985; Takeda 1991; Tomkin & Lambert 1983). Recently, Thevenin & 
Idiart (1999) studied NLTE effects in Fe abundances in metal-poor late-type stars. They found 
that ionized lines are not significantly affected by NLTE, and that NLTE corrections become less 
important as [Fe/H] increases, being minimal for solar abundance stars. Mostly high excitation 
lines are used in our analysis, so the NLTE effects that might be present should be minimal. 

Tomkin & Lambert (1983) found that NLTE effects cancel out when a differential analysis is 
carried out relative to a very similar star in terms of effective temperature and luminosity. Nearby 
stars of similar effective temperature and luminosity as the GC stars (see Figure 2) have been 
analyzed similarly in this paper, providing a differential abundance comparison that should remove 
any NLTE effects that might be present. 

4. Results 

Each iron abundance value for a particular Fe I line and star was derived from a comparison 
between the synthetic spectrum generated by MOOG and the observed spectrum. The [Fe/H] 
results for each line for cool, luminous stars in the solar neighborhood and GC stars are listed in 
Table 5 and 6, respectively. The uncertainties per line are estimated from the ability to distinguish 
models with different [Fe/H] considering the S/N of the observed spectrum. For cool, luminous 
stars in the solar neighborhood, the uncertainty per line is estimated to be ±0.05 dex, since the S/N 
is very homogeneous among all lines and stars. For GC stars, the uncertainty is listed individually 
for each line in Table 6. Figures 1, 7 and 8 show the synthetic and observed spectra for each line 
for a Ori, cool, luminous stars in the solar neighborhood, and GC stars, respectively. 

The final [Fe/H] results for each star are listed in Table 9. These values are the mean [Fe/H] 
weighted by the individual uncertainties in [Fe/H] for each Fe I line, as listed in Tables 5 and 6. 
The final uncertainty in Table 9 is the total uncertainty from Table 7 and 8. 

The analyses of our paper and Carr et al. (2000) have five stars in common: IRS 7, a Ori, HR 
6146 (30 Her), HR 6702, and HR 8383 (VV Cep). There is agreement in the Fe abundances within 
0.1 dex for four stars (IRS 7, a Ori, HR 6146, and HR 8383). Our Fe abundance for HR 6702 is 
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0.18 dex lower than that of Carr et al. (2000). The abundance technique used in both studies is 
very similar, but Carr et al. (2000) included a different set of Fe I lines and slightly different stellar 
parameters. 

IRS 7 is known to be a variable supergiant (Blum et al. 1996a). There is one Fe line common 
to the analysis of Carr et al. (2000) and of this paper, which was observed at different epochs 
by Carr et al. (2000) and by us. For this line, our result is consistent within the uncertainties 
with the results of Carr et al. (2000). Another spectrum of IRS 7, taken with NIRSPEC at the 
Keck Telescope one year after our observationes was kindly made available to us (Figer, private 
communication). We have one Fe line in common with the NISPEC spectrum, and again the iron 
abundance results from both data sets are also consistent within the uncertainties. We conclude 
that the effect of variability on our IRS 7 results are negligible considering our uncertainties. 

4.1. Galactic Center Mean [Fe/H] 

An unweighted mean value of [Fe/H] = +0.09 is obtained for the GC stars. A mean weighted 
by the total uncertainties in [Fe/H] (Table 9) gives a value of [Fe/H] = +0.12 for the GC stars. 
The estimated dispersion or a is 0.22 dex, which is not significantly different from the typical total 
uncertainty in [Fe/H] for each GC star of 0.28 dex. Even if the total uncertainties are overestimated 
by ~0.05 dex by the quadratic addition of the individual uncertainties, the typical uncertainty per 
star does not get small enough to resolve the [Fe/H] distribution (see Sec. 3.3). This means that 
the dispersion of the GC [Fe/H] distribution can be understood by the uncertainties present in the 
data. 

Note that the supergiant VR 5-7, at R = 30 pc, has [Fe/H] = +0.09 ± 0.22, which is similar 
to the supergiants at R < 0.5 pc, which have a mean [Fe/H] = +0.21 ± 0.15. 

An unweighted mean of [Fe/H] = -0.05 and a mean weighted by the total uncertainties of 
[Fe/H] = +0.03 was obtained for the eleven cool, luminous stars in the solar neighborhood. The 
estimated dispersion or a is 0.16 dex, which again is not significantly different from the typical 
total uncertainty in [Fe/H] which is 0.20 dex for these stars. 

The obtained mean [Fe/H] for the GC is close to solar, and furthermore, very similar to the 
cool, luminous stars in the solar neighborhood that were analyzed in the same way. 

5. Discussion 

Our mean [Fe/H] for the GC stars, +0.12 ± 0.22, is similar to that of supergiants in the 
Galactic disk. Figure 9a compares the GC [Fe/H] distribution with the [Fe/H] distribution of 40 
F, G, K, and M supergiants within 2 kpc of the Sun (Luck Sz Bond 1989). The mean [Fe/H] of the 
solar neighborhood supergiant stars from (Luck & Bond 1989) is +0.13 with a dispersion of 0.20 
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dex. The width and mean of the two distributions agree closely. In Figure 9b, we compare the 
GC [Fe/H] distribution with that of the 11 cool supergiants and luminous giants we have analyzed 
in this paper. Again, the mean and width of the [Fe/H] distribution for the GC stars and solar 
neighborhood stars are consistent. Unlike the Luck & Bond (1989) sample of supergiants, our 
11 nearby stars are restricted to the same range in stellar parameter space as the GC sample; in 
addition, we have used the same analysis method, Fe I lines, and stellar models for all the stars. It 
is important to emphasize that this differential analysis makes the abundance comparison between 
our samples of GC and nearby stars robust against possible systematic errors that might effect the 
absolute values of [Fe/H]. Hence, we can conclude that the Fe abundances for the sample of GC 
stars in this paper are nearly identical to those of similar T c g- and luminosity stars in the solar 
neighborhood, with no evidence for a super-solar metallicity. 

Our mean [Fe/H] for the GC stars of +0.12 ± 0.22 appears to be in conflict with the iron 
abundance of the Pistol star, at R = 30 pc in the Quintuplet cluster. Najarro et al. (1999) derive 
[Fe/H] ~ 0.5 for the Pistol star. They qualify their result as preliminary, since details of the 
atmosphere modeling, such as the effects of charge exchange reactions, have to be studied and 
included in future work. The Pistol star abundance is derived from emission lines in the winds 
from this hot star, so the atmospheric modeling techniques are fundamentally different from those 
more established techniques we have used to study abundances from photospheric absorption lines 
in cool giant or supergiant stars. 

Serabyn & Morris (1996) have proposed that the central cluster was built up by continuous, 
perhaps episodic, star formation over the lifetime of the Galaxy. In this case, the central 100 pc 
would form a stellar population distinct from the old population of the Bulge. In Figure 9c, the GC 
[Fe/H] distribution is compared with the [Fe/H] distribution for 262 stars in Baade's Window (BW) 
in the bulge (l,b = 1°, —4°; Sadler et al. 1996). Sadler et al. (1996) have a typical uncertainty per 
star of ±0.2 dex, so their finding that [Fe/H] ranges from —1.0 to +0.5 in BW shows the abundance 
spread of BW's distribution is well-resolved by their technique. The distribution of [Fe/H] in the 
GC, however, has a similar width as the uncertainties in [Fe/H], 0.28 dex per star. Stars as metal- 
rich or as metal-poor as observed in the Bulge are not found in our sample of GC stars. However, 
our GC stars are restricted to a young to intermediate age population, and hence an old stellar 
population, if it exists in the central 100 pc, is not sampled by the work of this paper. A larger 
abundance spread could be represented in such an older population, but a spectroscopic abundance 
analysis will require the sensitivity to reach K and M giants. A more complete comparison must also 
include measurements of the abundance patterns in the central regions relative to the Bulge. For 
example, the Bulge K giants analyzed by McWilliam & Rich (1994) show a distinct enhancement 
in the a— elements, compared to Fe. 

An IMF weighted toward more massive stars in the GC has been proposed by Morris & Serabyn 
(1996). Recently, Figer et al. (1999) have found an IMF slanted towards massive stars for the Arches 
and Quintuplet clusters, both located within 30 pc of the GC. A history of chemical evolution 
dominated by massive stars is expected to result in enhancements of oxygen and a— elements, such 
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as Mg, Si, Ca, and Ti, relative to Fe (Wheeler et al. 1989). The next step of our GC project is to 
measure abundances of a— elements and study selective enrichment in the GC. 

Chemical evolution models try to explain abundance patterns by considering the relative star 
formation rate, the gas infall and outflow rates, the star formation history and the abundance of 
the gas compared to the stars (see Audouze & Tinsley 1976; Rana 1991). Most chemical evolution 
models that reproduce the radial gradients in galaxies are valid only for distances greater than 2 
kpc from the Galactic Center (Tinsley & Larson 1978; Samland et al. 1997; Chiappini et al. 1999; 
Portinari & Chiosi 1999). Now that the molecular gas in the GC (from which stars form) has 
been extensively studied (Stark et al. 1991; Blitz et al. 1993; Morris & Serabyn 1996), and that the 
abundances of stars in the GC has been obtained for the first time, we strongly urge theoreticians to 
put all these data together into a detailed chemical evolution model describing the central parsecs 
of our Galaxy. 

6. Conclusions. 

We present the first measurement of stellar [Fe/H] for ten stars in the Galactic Center (GC), 
at distances from the GC of R < 30 pc. Nine GC stars are located in the central cluster (R < 0.5 
pc) and one star is located in the Quintuplet cluster (R = 30 pc). The abundance analysis is based 
on high-resolution (A/AA = 40,000) K— band spectra. The mean [Fe/H] of the GC is determined 
to be near solar, [Fe/H] = +0.12 ± 0.22, and also similar to the mean [Fe/H] for cool, luminous 
stars in the solar neighborhood, [Fe/H] = +0.03 ± 0.16, observed and analyzed in the same way. 
The width of the GC [Fe/H] distribution, which ranges from [Fe/H] = -0.3 to +0.5, is found to be 
significantly narrower than the width of the [Fe/H] distribution of Baade's Window, which ranges 
from [Fe/H] = -1.0 to +0.8. The GC [Fe/H] distribution is consistent with the [Fe/H] distribution 
of supergiant stars in the solar neighborhood. This suggest that the most luminous stars in the 
GC are unlikely to be dominated by bulge-like stars, and that the evolutionary path of the GC, 
while unique, is closer to the disk's than to the bulge's. The Quintuplet star at R = 30 pc has 
a similar [Fe/H] to stars located in the central cluster at R < 0.5 pc. In the future, abundance 
measurements of CNO and a— elements are planned to provide a complete view of the abundance 
patterns of the stars in the central regions of the Milky Way. 
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Table 1. Stellar Parameters for Nearby Late Type Stars. 



Star 


Spectral 


S/N a 


rp a 


M bol a 


log g a 








Type 




(K) 






(km s _1 ) 


(km s _1 ) 


HR 6146 


M6 III 1 


45 2 


3250±100 1 


-5.5 1 


O^iO.S 1 


2.0±0.5 2 


9.6±1.1 2 


HR 6702 


M5 II-III 1 


36 2 


3300IL100 1 


-3.4 1 


O^iO.S 1 


2.0±0.5 2 


7.8±1.8 2 


HR 7442 


M5 Illas 3 


36 2 


3450±100 3 


-4.0 3 


0.5±0.3 3 


2.4±0.2 2 


10.5±2.7 2 


HR 8062 


M4 Illas 3 


34 2 


3450±100 3 


-3.5 3 


0.7±0.3 3 


1.7±0.2 2 


8.9±0.5 2 


a Ori 


Ml-2 Iab-a 4 


305 5 


3540±260 4 


-7.4 6 


0.0±0.3 6 


2.8±0.2 2 


14.7±0.5 2 


HR 8383 


M2 Iape+ 7 


100 2 


3480±250 4 


-6.8 7 > 8 


0.0±0.3 7 


2.7±0.2 2 


14.4±2.1 2 


HD 202380 


M3 lb 7 


46 2 


3600±200 7 


-5.7 7 ' 8 


0.6±0.5 7 


2.5±0.2 2 


16.2±0.5 2 


HD 163428 


K5 II 7 


67 2 


3800±200 7 


-5.5 7 ' 8 


0.6±0.5 7 


2.3±0.2 2 


11.5±3.2 2 


BD+59 594 


Ml lb 9 


49 2 


4000±200 9 


-6.6 8 


-0.9±0.3 9 


3.1±0.2 2 


18.2±1.0 2 


HD 232766 


Ml lab 9 


38 2 


4000±200 9 


-6.6 8 


0.2±0.3 9 


2.2±0.2 2 


11.5±3.7 2 


HR 8726 


K5 lab 10 


57 2 


4000±200 10 


-5.2 10 > f 


! 0.5±0.3 10 


2.4±0.2 2 


11.4±0.5 2 



a S/N = mean signal to noise ratio per pixel ; T e ff = effective temperature; Mbd = bolo- 
metric magnitude; g = surface gravity; £ = microturbulent velocity; ( = macroturbulent 
velocity. 

References. (1) Smith & Lambert (1985); (2) this paper; (3) Smith & Lambert 
(1986); (4) Carr et al. (2000); (5) Wallace & Hinkle (1996); (6) Lambert et al. (1984); 
(7) Luck (1982b); (8) Lang (1991); (9) Luck & Bond (1989); (10) Luck (1982a). 
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Table 2. Fe Line Data. 



Wavelength (A) 


X (eV) a 


log gf a 


21781.82 


3.415 


-4.485 


22381.27 


5.844 


-1.458 


22386.90 


5.033 


-0.481 


22391.22 


5.320 


-1.600 


22398.98 


5.099 


-1.249 


22818.82 


5.792 


-1.296 


22838.60 


5.099 


-1.325 


22852.17 


5.828 


-0.612 



a X = Excitation Potential; gf = gf- 
value determined in this paper. 



-18- 



Table 3. Stellar Parameters for Galactic Center Stars. 



Star 


Luminosity 


S/N a 


T cff a 


M bol a > 6 


M a 


log g a 








Class 




(K) 




(M ) 




(km s _1 ) 


(km s _1 ) 


IRS 7 C 


I 


52 


3470±250 


-9.0 


17±3 


-0.6±0.2 


3.3±0.4 


20.6±2.7 


VR 5-7 


I 


30 


3500±300 


-7.8 


14±2 


-0.2±0.3 


2.9±0.5 


12.6±1.6 


IRS 19 


I 


75 


3650±300 


-7.2 


14±2 


0.1±0.3 


2.7±0.5 


13.4±2.0 


IRS 22 


I 


25 


3550±300 


-6.4 


10±2 


0.3±0.3 


2.5±0.5 


12.8±1.6 


BSD 124 


I 


10 


3600±300 


-5.5 


7±3 


0.4±0.3 


2.4±0.5 


12.7±2.7 


BSD 129 


I 


13 


3650±300 


-5.3 


7±3 


0.5±0.3 


2.4±0.5 


12.3±2.7 


BSD 72 


I 


13 


3750±300 


-4.5 


5±2 


0.8±0.3 


2.1±0.5 


11.0±2.7 


BSD 114 


III 


15 


3100±280 


-5.8 


3±1 


-0.2±0.5 


2.9±0.6 


9.0±3.0 


IRS 11 


III 


16 


3100±280 


-5.3 


3±1 


0.0±0.5 


2.8±0.6 


9.1±2.2 


BSD 140 


III 


13 


3100±280 


-4.8 


2.5±1 


-0.1±0.5 


2.9±0.6 


9.0±3.0 



a S/N = mean signal to noise ratio per pixel; T e g = effective temperature; M^oi = bolometric 
magnitude; M = mass; g = surface gravity; £ = microturbulent velocity; £ = macroturbulent 
velocity. 

b Error in M^ \ is ±0.4 (Carr et al. 2000), dominated by the uncertainty in the extinction 
curve (E(H -K)/A K ). 

c Stellar parameters from Carr et al. (2000). 
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Table 4. Fe I Equivalent Widths for Nearby Late- Type Stars. 









Wavelength of Fe I lines (A) 








22381.3 


22386.9 


22391.2 


22399.0 


22818.8 


22838.6 


22852.2 


Star 




Equivalent Widths of Each Fe I line (: 


mA) 




HR 6146 


90± 6 


355± 8 


162± 7 


210± 7 


130± 6 


277± 7 


211± 7 


HR 6702 


67± 5 


332± 9 


136± 7 


189± 7 


119± 6 


234± 7 


220± 7 


HR 7442 


60± 5 


339± 9 


143± 8 


174± 8 


Ul± 8 


294± 8 


246± 8 


HR 8062 


65± 5 


351± 9 


149± 7 


199± 6 


110± 6 


208± 7 




a Ori 


101± 8 


520±10 


218±10 


319±10 


178±10 


339±10 


271±10 


HD 202380 


102± 8 


476±11 


180± 9 


311± 9 


137± 8 


299± 9 




HR 8383 


86± 7 


466±10 


183± 9 


285± 9 


112± 8 


261± 9 


233± 9 


HD 163428 


87± 7 


430± 9 


172± 9 


241± 9 


113± 8 


263± 8 




BD+59 594 


107± 9 


522±11 


171±12 


329±12 


96± 9 


289±13 


226±12 


HD 232766 


66± 6 


410± 8 


159± 9 


228± 9 


97± 8 


216± 9 


173± 9 


HR 8726 


71± 6 


397± 8 


158± 8 


234± 8 


115± 8 


239± 8 


186± 8 
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Table 5. [Fe/H] for Nearby Late Type Stars. 



Wavelength of Fe I lines (A) 





21781.8 


22381.3 


22386.9 


22391.2 


22399.0 


22818.8 


22838.6 


22852.2 


Star 






[Fe/H] a Determined from Each Fe I Line 






HR 6146 


-0.05 


+0.05 


-0.20 


0.00 


-0.30 


+0.15 


+0.25 


+0.10 


HR 6702 


-0.05 


-0.05 


-0.25 


-0.05 


-0.30 


+0.20 


+0.10 


+0.30 


HR 7442 


+0.10 


+0.10 


-0.25 


+0.10 


-0.60 


+0.20 


+0.20 


+0.20 


HR 8062 


+0.10 


-0.10 


-0.05 


+0.05 


-0.20 


+0.10 


-0.05 




a Ori 


+0.15 


0.00 


+0.05 


+0.05 


-0.10 


+0.20 


+0.10 


0.00 


HD 202380 


+0.10 


+0.10 


+0.10 


0.00 


+0.10 


+0.10 


+0.10 




HR 8383 


-0.05 


-0.10 


-0.10 


-0.10 


-0.20 


-0.15 


-0.25 


-0.10 


HD 163428 


+0.07 


0.00 


0.00 


-0.02 


-0.25 


-0.05 


-0.05 




BD+59 594 


-0.10 


0.00 


-0.25 


-0.25 


-0.25 


-0.30 


-0.35 


-0.35 


HD 232766 


-0.10 


-0.20 


-0.15 


-0.15 


-0.40 


-0.20 


-0.40 


-0.40 


HR 8726 


+0.15 


-0.10 


0.00 


-0.05 


-0.20 


0.00 


-0.10 


-0.20 



a Error in [Fe/H] for each line is ± 0.05 dex. 
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Table 6. [Fe/H] for Galactic Center Stars. 



Wavelength of Fe I lines (A) 





21781.8 


22381.3 


22386.9 


22391.2 


22399.0 


22818.8 


22838.6 


22852.2 


Star 






[Fe/H] Determined from Each Fe I Line 






IRS 7 


+0.35 


±0.22 


— U.lO 


— U.1U 




+U.oo 


-0.30 


— 




±0.05 


±0.07 


±0.10 


±0.05 




±0.05 


±0.05 




VR 5-7 

V 11 U 1 


-1-0 1 2 


-1-0 20 


-0.30 


±0.15 


-0.40 


±0.10 


-1-0 1 5 






±0.03 


±0.07 


±0.10 


±0.05 


±0.10 


±0.05 


±0.07 




IRS 19 


±0.30 


±0.40 


-0.30 


±0.40 


-0.05 


±0.45 


-0.05 






±0.05 


±0.05 


±0.10 


±0.07 


±0.10 


±0.05 


±0.10 




IRS 22 


±0.30 


-0.05 


-0.40 


±0.15 


-0.40 


±0.20 


±0.35 


±0.30 




±0.05 


±0.05 


±0.10 


±0.10 


±0.10 


±0.10 


±0.10 


±0.10 


BSD 124 




±0.22 


±0.20 


±0.20 


-0.10 


±0.20 


±0.20 








±0.07 


±0.10 


±0.20 


±0.10 


±0.05 


±0.10 




BSD 129 




±0.50 


±0.50 


±0.70 


±0.20 


±0.60 


±0.60 








±0.20 


±0.10 


±0.20 


±0.10 


±0.10 


±0.10 




BSD 72 




±0.10 


±0.40 


±0.35 


±0.05 


±0.35 


-0.20 








±0.10 


±0.10 


±0.10 


±0.10 


±0.10 


±0.10 




BSD 114 




-0.20 


-0.85 


-0.10 


-0.70 


0.00 


-0.30 








±0.10 


±0.20 


±0.10 


±0.10 


±0.10 


±0.10 




IRS 11 


-0.15 


±0.25 


-0.80 


-0.20 


-0.80 


-0.15 


-0.50 






±0.05 


±0.05 


±0.10 


±0.10 


±0.10 


±0.10 


±0.10 




BSD 140 




0.00 


-0.40 


±0.10 


-0.40 


0.00 


0.00 








±0.20 


±0.20 


±0.20 


±0.20 


±0.10 


±0.10 
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Table 7. Uncertainties in [Fe/H] for Solar Neighborhood Stars. 



Star 


±T cfr a 


±log g a 


±e a 


±c a 


Std. b 


Total c 


HR 6146 


T0.03 


±0.10 


to.h 


±0.08 


±0.07 


±0.19 


HR 6702 


T0.03 


±0.10 


to.h 


±0.13 


±0.08 


±0.22 


HR 7442 


T0.03 


±0.10 


T0.04 


±0.19 


±0.11 


±0.25 


HR 8062 


T0.03 


±0.10 


T0.04 


±0.04 


±0.04 


±0.13 


a Ori 


T0.07 


±0.10 


T0.04 


±0.04 


±0.03 


±0.14 


HD 202380 


T0.05 


±0.17 


T0.04 


±0.04 


±0.02 


±0.19 


HR 8383 


T0.07 


±0.10 


T0.04 


±0.15 


±0.03 


±0.20 


HD 163428 


T0.05 


±0.17 


T0.04 


±0.22 


±0.04 


±0.29 


BD+59 594 


T0.05 


±0.10 


T0.04 


±0.07 


±0.05 


±0.15 


HD 232766 


T0.05 


±0.10 


T0.04 


±0.26 


±0.05 


±0.29 


HR 8726 


T0.05 


±0.10 


T0.04 


±0.04 


±0.04 


±0.13 



a T c g = effective temperature; g = surface gravity; £ = microtur- 
bulent velocity; ( = macroturbulent velocity. 

b standard error, determined from scatter among [Fe/H] values 
measured for different Fe I lines. 

c total uncertainty, derived from quadratic sum of uncertainties 
due to the standard error and the uncertainties from the stellar 
parameters (T eS , log g, £, (). 
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Table 8. Uncertainties in [Fe/H] for Galactic Center Stars. 



Star 


±T cfr a 


±log g a 


±f a 


±c a 


Std. b 


Total c 


IRS 7 


T0.07 


±0.07 


T0.09 


±0.19 


±0.13 


±0.27 


VR 5-7 


T0.08 


±0.10 


to.h 


±0.11 


±0.11 


±0.23 


IRS 19 


T0.08 


±0.10 


to.h 


±0.14 


±0.13 


±0.26 


IRS 22 


T0.08 


±0.10 


to.h 


±0.11 


±0.12 


±0.24 


BSD 124 


T0.08 


±0.10 


to.h 


±0.19 


±0.06 


±0.26 


BSD 129 


T0.08 


±0.10 


to.h 


±0.19 


±0.08 


±0.27 


BSD 72 


T0.08 


±0.10 


to.h 


±0.19 


±0.10 


±0.27 


BSD 114 


T0.08 


±0.17 


T0.13 


±0.21 


±0.16 


±0.35 


IRS 11 


T0.08 


±0.17 


T0.13 


±0.15 


±0.18 


±0.33 


BSD 140 


T0.08 


±0.17 


T0.13 


±0.21 


±0.10 


±0.32 



a T c ff = effective temperature; g = surface gravity; £ = micro- 
turbulent velocity; £ = macroturbulent velocity. 

b standard error, determined from scatter among [Fe/H] values 
measured for different Fe I lines. 

c total uncertainty, derived from quadratic sum of uncertain- 
ties due to the standard error and the uncertainties from the 
stellar parameters (T c g, log g, £, £)• 
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Table 9. Mean [Fe/H] for each star. 



Solar Neij 


^hborhood Stars 


Galact 


ic Center Stars 


Star 


Nlines 


[Fe/H] a 


Star 


Nlines 


[Fe/H] a 


HR 6146 


8 


-0.01±0.19 


IRS 7 


6 


+0.09±0.27 


HR 6702 


8 


-0.02±0.22 


VR 5-7 b 


7 


+0.09±0.23 


HR 7442 


8 


+0.01±0.25 


IRS 19 


7 


+0.29±0.26 


HR 8062 


7 


-0.03±0.13 


IRS 22 


8 


+0.09±0.24 


a Ori 


8 


+0.05±0.14 


BSD 124 


6 


+0.17±0.26 


HD 202380 


7 


+0.07±0.19 


BSD 129 


6 


+0.49±0.27 


HR 8383 


8 


-0.14±0.20 


BSD 72 


6 


+0.17±0.27 


HD 163428 


7 


-0.05±0.29 


BSD 114 


6 


-0.29±0.35 


BD+59 594 


8 


-0.24±0.15 


IRS 11 


7 


-0.16±0.33 


HD 232766 


8 


-0.26±0.29 


BSD 140 


6 


-0.06±0.32 


HR 8726 


8 


-0.07±0.13 









a Error in [Fe/H] is the total uncertainty from Table 7 and 8. 
b M supergiant star in the Quintuplet cluster (R = 30 pc). 
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Fig. 1. — Observed spectrum (filled squares) of a Ori from Wallace & Hinkle (1996) compared with 
synthetic spectra generated by MOOG, using CNO abundances from Lambert et al. (1984), model 
atmosphere from Plez (1992), and stellar parameters listed in Table 1. The error bars come from 
the difference of observations taken in two epochs. We overplot a synthetic spectrum (thick line) 
derived using the Fe abundances listed in Table 5 and synthetic spectra (thin lines) computed 
with Fe abundances different by ±0.20 dex , which corresponds to the typical uncertainty in [Fe/H] 
for a solar neighborhood cool, luminous star. Fe lines and their wavelengths in A are marked 
(bold vertical lines) at the top of each panel. CN, Sc, and unidentified lines (question marks) are 
also marked (vertical lines). Tickmarks along the x-axis are 1 A apart. The [Fe/H] value from 
Table 5 for each line is given in each panel. 

Fig. 2. — CO index vs effective temperature, T e g , for supergiant and giant stars. Upper panel: the 
calibration for supergiant stars, with T e g from Dyck et al. (1998) (filled squares) and Richichi 
et al. (1998) (open circles), both measured from lunar occultations. Lower panel: the calibration 
for giant stars. Values of T c g from McWilliam (1990) (open triangles), Smith & Lambert (1985, 
1986, 1990) (filled circles), and Fernandez- Villacahas et al. (1990) (open stars) are derived from 
a relationship between V — K and T e g measured by lunar occultation techniques. T c g from Dyck 
et al. (1998) (open circles), are derived from lunar occultation measurements. In both panels the 
uncertainties come from the literature. The best unweighted linear fit to the data (solid lines) is 
used to compute the T e g- of the galactic center stars. 

Fig. 3. — HR diagram of galactic center (GC) stars (filled circles) and cool, luminous stars in the 
solar neighborhood (open squares). Stellar parameters and their uncertainties are given in Tables 
1 and 3. Solar metallicity evolutionary tracks (Schaller et al. 1992) are overplotted to determine 
masses for GC stars. The GC stars and the cool, luminous stars in the solar neighborhood occupy 
the same place in the HR diagram, hence are very similar types of stars. 

Fig. 4. — Relation between microturbulent velocity (£) and surface gravity (log g). Values of £ 
and log g are from Table 1 of this paper (open squares), Carr et al. (2000) (crosses), and Smith 
& Lambert (1985, 1986, 1990) (filled triangles). The best unweighted linear fit to the data 
(solid line) is used to compute £ for the GC stars. 

Fig. 5. — Left panels: \ 2 contour maps for determining the nitrogen abundance, log e (N), and 
macroturbulent velocity (() for HR 8726. The minimum x 2 gives the best set of log e (N) and £ 
for each set of CN lines (21798 A - 21804 A at the top, and 22403 A - 22409 A at the bottom). 
Right panels: observed spectrum of HR 8726 (filled squares) with error bars from the difference of 
observations taken at two positions along the slit and the synthetic spectrum of HR 8726 (solid line) 
using the derived log e (N) and ( for each CN band (21798 A - 21804 A at the top, and 22403 A - 
22409 A at the bottom). The synthetic spectrum is calculated by MOOG, using stellar parameters 
in Table 1 and a model atmosphere from Plez (1992). The value of log e (N) derived by this 
technique is not a true nitrogen abundance, because a combination of C and N abundances are 
needed to derive a true value of log e (N) from CN lines. 
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Fig. 6. — Relation between macroturbulent velocity (Q and surface gravity (log g) for supergiant 
stars {filled circles). All supergiant stars from Table 1 and IRS 7, VR 5-7, IRS 19, and IRS 22 
from Table 3 are plotted. The best unweighted linear fit to the data {solid line) is used to compute 
C for the GC supergiant stars BSD 72, BSD 124, and BSD 129. 

Fig. 7. — Observed {filled squares) and synthetic spectra {solid line) for cool, luminous stars in 
the solar neighborhood. The synthetic spectrum computed by MOOG, using the stellar parameters 
listed in Table 1 and model atmospheres from Plez (1992). The error bars come from the difference 
of observations taken at two positions along the slit. Fe lines and their wavelengths in A are marked 
{bold vertical lines) at the top panels. CN, Sc, and unidentified lines {question marks) are also 
marked {vertical lines). Tickmarks along the x-axis are 1 A apart. The [Fe/H] value from Table 
5 for each line is given in each panel. 

Fig. 8. — Similar to Fig. 7, except these are observed and synthetic spectra for galactic center 
stars, using stellar parameters given in Table 3. Symbols are the same as Fig. 7. The [Fe/H] value 
from Table 6 for each line is given in each panel. 

Fig. 9. — (a) Fractional distribution of [Fe/H] for 10 GC stars ( solid line) compared to [Fe/H] 
for 40 solar neighborhood supergiant stars {dashed line) from Luck & Bond (1989). (b) Fractional 
distribution of [Fe/H] for 10 GC stars {solid line) compared to [Fe/H] for 11 solar neighborhood 
stars observed and analyzed in this paper {dashed line), (c) Fractional distribution of [Fe/H] for 10 
GC stars {solid line) compared to [Fe/H] for 262 Baade's window stars {dashed line) from Sadler 
et al. (1996). 



- 31 - 



- ^ (\ 


1 


T 1 O- "V\ 








CO 




- C\? y 




CV 


O 


" CD / 


O 




+ _ 







o 

CO 

cv 

C\2 









- M . 






c\i 






- lO 


t>- ^^?4L. 




00 






- C\2 






CV 




O 






O - 


CD 




o 








iff , 1 


1 





CO 

o 



CO 

o 



- o 


1 

o 




1 




CO 






CT) 








05 








CO 








" CV 






o 


CV 






T 1 


CD 


^ 1 

1 \ 




o 


- 


1 


1 

1 




i — 1 


CO 


CO 






o 


o 



Y 



32 



4000 



3000 



CD 



E ' 



4000 



3000 



i 1 1 r 



i r 



■ Dyck 

O Richichi 



A 



Supergiant stars 



Giant stars 



A McWilliam 

■ S&L 

- -A Fernandez 

O Dyck 




J I I L_ 



J L_ 



10 



20 



CO index 



- 33 - 




-34- 



i r 



i i r 



i i r 



i i i r 



□ this paper 

x Carr et ah 
A Smith & harnbert 







4 1 




i ^ 


\^ i- 


> E 


3 




GC range 



J I I I I I I I I I I I I I I I I I I I I I I I I I L 

-1 -0.5 0.5 1 

log g 



, I , L_ 

8 10 12 

<f(km/s) 



0.85 



t= 11.0 



14 



21800 



21804 



A (A) 



22403A-E2-4-O-&A— 

IIR 8726 




X5 

o 

N 



S-H 

o 



0.9 



0.85 



8 



10 12 

((km/ s) 



14 



1 


1 


CN CN CN*" , CN CN 


CN CN/ 














• 






1 


10.6 
i 1 


22404 

A (A) 


22408 



- 36 - 




-37- 




h CO CD h CO tD CO CD CO CD 00 CD 



oo oo oo oo oo 



sqoj pszijieiiiJON 



-38- 




S qoJ pszijieiiiJON 



-39- 



z' 
u 





C/J 



"El," 



Z_ 

u 



cu 

"fen 



C0- 

o 

I, 



o 
d 





CO 

o 



CD 
O 



CO 

O 



CO 

o 



CO 

o 



CD 
O 




T 

H co 
cv 
cv 




o- 
d 




co 



cv 

Q 

m 
CO 



_L- 



CO CD 
O O 



CO CD 

O O 



sqoj pszijieiiiJON 



COCO h [Q tD CO CD ^HCOCD CO CD 

OO OO OO OO OO 

S qoJ pszijieiiiJON 



41 



0.6 



0.3 







0.6 



o 



^0.3 




0.6 



0.3 







GC (N Total = 10) 
SG (N Total = 40) 



gc (: 

SN (N 



Total 



Total 



10) 

11 



GC (N 
BW (N 



Total 



'Total 



10) 
262) 



(a) 



(b 



(c 



2 -1.6 -1.2 -0.8 -0.4 0.4 0.8 1.2 

[Fe/H] 



